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We investigate the phase diagram in the plane of temperature and chemical potential 
mismatch for an asymmetric fermion superfiuid with double- and single-species pairings. 
There is no mixing of these two types of pairings at fixed chemical potential, but the 
introduction of the single species pairing cures the magnetic instability at low tempera- 
ture. 



1. Introduction 

Inspired by ultracold atomic physics, nuclear physics and color superconductiv- 
ity, the fermion pairin g between different species with mismatched Fermi surfaces 
prompted great intereslP^ ^ * 4 * 5 * 6 * 7 * 8 * 9 * 10 * 11 * 1 ' 2 * 13 ! in recent experimental and theo- 
retical studies. In conventional fermion superfiuid the ground state is well described 
by the BCS theory, while for asymmetric fermion superfiuid with double-species 
pairing the phase structure is much more rich and the pairing mechanism is not 
yet very clear. Var ious exot ic phases have been suggested in the literatures, such 
as the Sarma 

phasJSGSEEl the Fulde-Fe rrel-L arkin-Ovchinnikov(FFLO) phase^ , 
the phase with deformed Fermi surfacefPESl 

and the phase separatio n I I in 

coordinate space. 

For asymmetric fermion superfiuid, one of the most surprising phenomena is the 
magnetic instability, i.e., when the asymmetry characterized by the chemical po- 
tential mismatch is larger than the condensate of Cooper pair, the BCS superfiuid 
is unstable under the perturbation of an external magnetic field. This instability 
leads to the ground state to be a FFLO phase. For color superconductivity, a sim- 
ilar but more complicated phenomenon is also found, namely the chromomagnetic 
instability^! One of the mechanisms to cure the magnetic instability is to introduce 
another pairing channel. In this paper, we propose a simple model including both 
double- and single-species pairings and investigate how the single-species pairing 
influences the FFLO state. 

Our model can be used to describe a general fermion superfiuid system where 
fermions from the same species can form Cooper pairs. In ultracold atom gas like 6 Li 
and 40 K, there exist pai rings between different elements and between different st ates 
of the same elemenlP2l. In color superconductivity phase of dense quark matter^l, 
the quarks of different flavors can form total spin zero Cooper pairs and the quarks 



1 



2, 2008 9:25 WSPC/INSTRUCTION FILE QM06-huang 



2 Xuguang Huang, Xuewen Hao and Pengfei Zhuang 

of the same flavor can combine into total spin one pairs too. In neutron stars 
proton-proton, neutron-neutron and neutron-proton pa irs are all po ssible . Recently 
discovered high temperature superconductivity of MgBjP^SEZESl can a j so we ^ 
described by an extended two-band BCS theory where the electrons from the same 
energy bands form Cooper pairs. 

2. The Model 

We consider a fermion system containing two species a and b with masses m a 
and mi, and chemical potentials fi a and fib satisfying fi a < fib- The system can be 
described by the Hamiltonian 



P,k L 



o P -P -k a k + 5a a p a -p a -k a k + 9bO p O_ p 0_ k k 



(1) 



where p and k are fermion momenta, a, b, a + and b + arc the annihilation and 
creation operators, the coupling constants g,g a and gi, are positive to keep the 
interactions attractive, e a — p 2 /(2m a ) — fi a and e/, = p 2 /(2mb) — fib are the 
non-relativistic dispersion relations, and in continuous limit we simply replace 
J2 P w hh J d 3 p/(2n) 3 . The (pseudo-)spin cr, a' =T, j has been introduced to keep 
our model satisfying the Pauli rule explicitly. Other inner structures like isospin, 
flavor and color are neglected since they are not essential for our purpose. The 
Hamiltonian has the symmetry U a (l) £§> Ub{l) with the element U(9 a , 9b) defined as 
U(9 a ,9 b )a p U+{9 a ,9 b ) = e w °a p and U(9 ai 9 b )b p U + {9 a ,9 b ) = e l6 4 p . 

In the framework of mean field approximation, after taking a Bogoliubov-Valatin 
transformation from fermions a and b to quasi-particles A and B, the thermody- 
namic potential can be expressed in terms of the quasi-particles, 

n = Al/g a + A 2 b /g b + 2A 2 /g 

+ ]T [e a + e b E A -E B - 2Tln (((1 + e ^/ r ) (l + e- E -' T ))] , (2) 
p 

where 

A a =g a J2(^-A), A 6 = 56 £<&L p S],>, A -"/ 2 E(' : » ( 3 ) 

P P P<? 

are the corresponding a-a, b-b and a-b pairing condensates, and 



4 = 4 + vW+(4f. E% = el-J(el) 2 + (el? (4) 
are the quasi-particle energies with 

4 = [(4 + A 2 +A 2 )±( e 2 + A 2 + A 2 )]/2, 
(el) 2 = A 2 \(e a - e b f + (A a - A b f} . (5) 
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A a =/= and A& ^ correspond, respectively, to the spontaneous symmetry breaking 
patterns t/ (l) £g) £4(1) — > £4(1) an d £4(1), and A ^ means the breaking pattern 
U a (l) <8> £4(1) -> £4-&(l) with the element £/(0) defined as U(9)a p U+(9) = e i9 a p 
and U(9)b p U + (9) = e~ l9 b p . For simplicity, we have considered all the condensates 
as real numbers. 

To study the FFLO state, we modify the definitions of the fermion dispersions 

as 

ef = (p± qi) 2 /(2mj) - m, i = a,b (6) 

with + corresponding to spin-up and — to spin-down, where are the FFLO 
momenta of condensates of z — z pairing, which together with the condensates are 
self-consistently determined by the coupled set of gap equations 

dSl/dA a = 0, dn/dA b = 0, dtt/dA = 0, 0fi/0<fc = (7) 

and the minimum thermodynamic potential. Note that the FFLO state we consid- 
ered here is the simplest form where the Cooper pairs in coordinate space have the 
plane wave forms 

A(x) = Ae^+^)- x , A Q (x) = A a e 2iq «' x , A b (x) = A b e 2 "^ x . (8) 

We will choose q a = q b = q, such a choice can avoid possible dyna mic insta bility 
due to the different superflow velocities of a — a and 6—6 channel s 11 Ob- 
viously, the translational symmetry and rotational symmetry in the FFLO state 
are spontaneously broken. To have a simple analytical expression for the spectral 
of quasi-particles, we set m a = m& = rn. In this case, the thermodynamic potential 
reads, 

n = A 2 Jg a + A 2 b /g b + 2A 2 /g + ]T [c~ + e 6 ~ - (E+ + E A + E+ + E B ) /2 

p 

-Tln((l + e- B l/ T ) (l + e^/ T ) (l + e^/ T ) (l + e^/^)], ( 9 ) 

where 



E% = \/e% + 6e* - + (e\r ± Se (10) 

are the quasi-particle energies with the notations 

4 = [(etc- +AI + A 2 ) + (e+e 6 " + A 2 + A 2 )] /2, 
4 = [(e+e- + A 2 + A 2 ) - (e+e 6 - + A 2 + A 2 )] /2, 
Se=(et-e-)/2=(e+-e b -)/2, 
(elf = A 2 [(e+ - e+)(e- c >) + (A 6 - A a ) 2 ] . (11) 
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The magnetic stability, namely the stability of the homogeneous superfluid 
against the pair momentum fluctuations is characterized by the superfluid density 

Ps = d 2 n/dci\ =0 . (12) 

Positive p s means stable homogeneous superfluid and negative p s indicates possible 
FFLO state. In this sense q plays the role of an external magnetic potential eA. 

In the end of this section, we calculate the gapless nodes for the case with zero 
FFLO momentum where at least one of or E% crosses the momentum axis, 

E+E A E+E B = (A a e b - A b e a f + (A 2 + A a A b + e a e b ) 2 = 0. (13) 

There are two classes of gapless solutions. 1) Only A is nonzero. The gapless mode 
happens at momenta 

p± - ^2m(//± v^ 2 -A 2 ), (14) 

where we have introduced the average chemical potential fx = (p a + p b ) /2 and the 
chemical potential mismatch Sp = (p b — p a ) /2. It is clear that the condition to 
have the gapless mode is Sfi > A. 2) All the three condensates are zero or only one 
of A a and A;, is nonzero. In these cases, we have e a = or e b = or e a e b — which 
mean real fermion excitations exactly at the Fermi surfaces. 

3. Phase Diagrams 

Since the model is non-renormalizable, we introduce in the numerical calculations 
a cutoff A = Xpp with A = 2 1//4 and the average Fermi momentum pp = \J2rap. 
We choose ppa — 0.4, ppa a = PFO-b — 0.32 with a — mg/(Air), a a = mg a /{AiT) and 
a b = mg b /(AiT) being the s-wave scattering lengths. We have checked that in the 
parameter region of < pf&, p_fo q , PF^b < 1 and < a a — a b < a there is no 
qualitative change in the phase diagrams. 

The phase diagram in T — 5p plane is shown in FigfT] The upper panel is in 
the familiar case without single-species pairing . The homogeneous BCS state can 
exist at low temperature and low mismatch, and the inhomogeneous FFLO state 
survives only in a narrow mismatch window. When the single-species pairing is 
included as well, see the bottom panel of Fig[T] the inhomogeneous FFLO state of 
a-b pairing is eaten up by the homogeneous superfluid of a-a and b-b pairings at 
low temperature, just as we expected, and survives only in a small triangle at high 
temperature. The phase transition from the a-b pairing superfluid to the a-a and 
b-b pairing superfluid is of first order. Note that for systems with fixed chemical 
potentials there is no mixed phase of double- and single-species pairings, and the 
situation is similar to a three-component fermion system^. 

4. Summary 

We have investigated the phase structure of an asymmetric two-species fermion 
superfluid with both double- and single-species pairings. Since the attractive inter- 
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Fig. 1. (color online) The phase diagram in T — <5/x plane at fixed fi = 50Ao with Ao being the 
symmetric gap at S^i = 0. The upper panel is with only double-species pairing, and in the bottom 
panel the single-species pairing is included as well. 



action for the single-species pairing is relatively weaker, it changes the conventional 
phase diagram at low temperature. In any system with chemical potential imbal- 
ance, the double-species pairing in FFLO state is replaced by the single-species 
pairing. In the region with single-species pairing, the interesting gaplcss superfluid 
is washed out and all fermion excitations are fully gapped. We should note that in 
this paper we considered only the grand canonical ensemble. For the case of canon- 

lOOl 

ical ensemble, the phase diagram becomes much more riclr^ 1 . 
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